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We find a first order Coulomb-Higgs phase transition at moderately large values of the coupling A, and no 
evidence for a change of order at any finite value of it. 
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1. INTRODUCTION 

The Higgs mechanism is an essential part of the 
present day formulation of the Standard Model. 
The U(l) gauge model coupled to scalars is a sim- 
plification of the gauge-Higgs sector of the SM, 
which yet keeps the unsolved problem of defin- 
ing non perturbatively a non-asymptotically free 
field theory, and which can be addressed with 
rather simple technical means. 

2. THE COULOMB HIGGS PHASE 
TRANSITION IN THE /3 k MODEL 

The three parameter U(l)-Higgs model is de- 
scribed by the action 
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In the A — cxD limit, 
simplifies to 



1 and the action 
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The phase diagram of that restricted version [1] 
has been discussed by us in Latt93 [2], and the dif- 
ficulty to identify the order of the Coulomb-Higgs 
phase transition was made apparent. The [3 ^ oo 
end of that line is the X-Y model phase transi- 
tion, a well established second order one, which 
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has induced researchers to think the line to be a 
prolongation of its end point, and as such, second 
order [4]. Yet, the relationship of the model to 
the q > 1 version and to the Zjv model reported 
in [2] and [3] points towards a possible first order 
transition. 

3. THE THREE PARAMETER MODEL 

The difficulty to study directly the A — cxd limit 
has lead us to study the problem at finite, variable 
A. This had been done earlier [5], and we intend 
to improve on the statistics. In order to achieve 
that goal, we have fixed [3 = 1.15, which is on the 
Coulomb-Higgs side of the first order Confining- 
Coulomb phase transition. We have worked at 
values of A 0.003, 0.005, 0.01, 0.03, 0.1, 0.3 on lat- 
tices 6'*,8'*,12'* and IG'*, with statistics ranging 
between 10® and 10® Montecarlo iterations per 
measurement, on workstations, on parallel ma- 
chines at the Edinburgh Parallel Computing Cen- 
tre and on our custom, 64 T800 processor com- 
puter RTN [6]. We have used an over-relaxed 
Metropolis update method, with a multicanoni- 
cal procedure to accelerate flip-flop rate when the 
latter is too low. The jack-knife method has been 
used in the error estimation. 

4. THE ENERGY HISTOGRAMS 

Fig. 1 shows typical histograms of the link en- 
ergy ($[/$), whose most striking features are: 

• Two-peak structure. 

• Assymetry of the peaks, with a narrow low- 
energy peak and a broad high-energy one. 
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• Strong size dependence, consisting on a nar- 
rowing of the gap and width of the peaks as 
the lattice size L increases, effected by a 
desplacement towards lower energy of both 
peaks, especially the broad one. 



histogram has been shifted to the apparent crit- 
ical point by the spectral density method [7]. A 
cubic spline has been preferred to other functional 
shapes because of its ability to reproduce a max- 
imum and accomodate the mixed states, whose 
influence on the histogram is otherwise difficult 
to account for. 
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Figure 1. Histograms of the link energy in an S'* 
lattice and a 12'* lattice at A = 0.01 



5. MEASUREMENTS OF THE LATENT 
HEAT 

In order to quantify the general features listed 
in the previous section, we have measured the la- 
tent heat at the apparent critical point at a given 
L, defined as the value Kc for which the link- 
energy spread (specific heat) is maximum. The 
definition is global, and so, little dependent on lo- 
cal details, and has a clear thermodynamic limit. 

The latent heat has been defined as the differ- 
ence between the positions of the maxima for a fit 
to each peak separately to a cubic spline, after the 
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Figure 2. Latent heat as a function of l/V at 
different values of A 



Fig. 2 shows, for A < 0.03, evidence for a non- 
zero-point extrapolation of the latent heat in the 
V ^ CO limit, suggesting the first order of the 
phase transition for smaller values of A. For A > 
0.03 the figure is not conclusive, so that on the 
basis of the behaviour of the latent heat alone, 
the ocurrence of a tricritical point at A > 0.03 
cannot be descarded. 

6. THE MAGNETIC SUSCEPTIBILITY 

In analogy with magnetic systems, we define 
the magnetic susceptibility as x = c^V , being 
the variance of each peak in the histogram, a is 
measured by fitting the log of the energy around 
each peak to a parabola, and the dependence of 
.^X on l/V is shown in Fig. 3. 

The figure shows a rather fiat behaviour of the 
susceptibility as — cxd at all values of A, with no 
trace of divergence as A increases. This behaviour 
seems to exclude the existence of a second order 
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Figure 3. Square root of the susceptibility of the 
higher state versus l/V at different values of A 



phase transition in any limit, pointing to the pos- 
sibility of the first order character being kept in 
the thermodynamic limit as A — cxd. 

7. FINITE SIZE EFFECTS 

The finite size effects observed can be ac- 
counted for almost completly by the shift of the 
apparent critical point produced by the change of 
lattice size. Fig. 4 shows the position of the high 
energy peak when the critical point, Kc, is deter- 
mined as the value for which the specific heat is 
maximum, which are the values used all through 
the paper, or as the value at which the two peaks 
have equal hight, another popular choice for the 
definition of the apparent critical point. Also 
plotted is the position of the peak when, at dif- 
ferent lattice sizes L, the value of Kc measured at 
_L = 16 is used. As can be seen, the position of 
the peak does not change practically. 

8. CONCLUSIONS 

The analysis of the latent heat shows that the 
Coulomb-Higgs phase transition is first order for 
all the finite values of A studied. The added ev- 
idence of the magnetic susceptibility not grow- 
ing with the lattice size and with A seems to ex- 
clude the existence of a finite Ac beyond which 
the transition becomes second order and points 
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Figure 4. Positions of the high energy peak as 
three different definitions of Kc (maximum specific 
heat, equal hight and Kc at L = 16 ) are used. A 
is 0.01 



towards the possibility of it being first order even 
at A = cxD (fixed modulus case) . 
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